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Abstract 

We study the properties of a few-electron system confined in coupled elongated quantum dots 
(QDs) using a model Gaussian potential and the numerical exact diagonalization technique. In 
the absence of magnetic fields, as the aspect ratio r between the QD extensions in the direction 
perpendicular and parallel to the coupling directions increases, the exchange energy exhibits a sharp 
variation at the specific value r = 3.9, before (after) which the exchange energy increases (declines). 
The sharp variation occurs because of a sudden change in the single-particle configuration of the 
triplet state. The stability region with one electron in each of the QDs is found to shrink, and 
finally vanishes as it becomes progressively easier to localize both electrons into the lower QD. For 
r > 3.9, the first singlet-triplet transition shifts to a small magnetic fields. 

PACS numbers: 73.21.La, 73.21.-b 
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I. INTRODUCTION 



Coupled quantum dots (QDs) based on two-dimensional electron gas (2DEG) formed 
with GaAs/AlGaAs heterostructures are promising candidates for quantum logic applica- 
tions because of the ability to coherently manipulate the many-body spin states by using 
external electromagnetic fields.™*^ Recently, a coherent controlled cycle of many-body state 
preparation, spin-interaction and projective read-out has been achieved in laterally coupled 
QDs. 4 In such an experiment, the electromagnetic control of the exchange energy J, which 
drives the Rabi oscillations between the lowest singlet and triplet states, is of utmost impor- 
tance. It is well known that the hyperfine interaction between the electron and nuclear spins 
competes with the exchange energy to destroy the singlet-triplet coherence.- Therefore, in 
order to retain the spin-state coherence in coupled GaAs/AlGaAs QDs, it is important to 
optimize the exchange energy to exceed the hyperfine interaction significantly. 

A wealth of theoretical work has been devoted to study the exchange energy in coupled 
QD systems.- 1 ^ The main focus of these studies is the tunability of the exchange energy by the 
electromagnetic fields and/or the parameters defining the interdot coupling strength, e.g., 
interdot separation and barrier height. The optimization of the exchange coupling J — given 
a fixed interdot distance, which is predetermined by the lithography of the top gates — has 
been rarely discussed. In this work, we investigate such a possibility by considering QDs 
elongated perpendicularly to the coupling direction. In this configuration, one can expect the 
overlap between the electron wavefunctions in the two QDs to increase, which will enhance 
their interactions. Our work is encouraged by the recent proposal of using coupled elongated 
QDs to construct robust spin-qubits with all-electrical qubit manipulation capabilities. 8 

In this paper, we perform a detailed analysis of the two-electron system in coupled elon- 
gated QDs to show that the exchange coupling indeed becomes larger with increasing aspect 
ratio between the extensions of each QD perpendicular and parallel to the coupling direc- 
tion (r = Ry/R x ). Our analysis based on the numerical exact diagonalization technique 
indicates that the cause of this enhancement is far from intuitive, while there is an optimum 
r value beyond which the exchange energy J decreases. Furthermore, for r > 5, we find 
that the stability region for one electron in each QD shrink to vanish. Finally, the magnetic 
field, which defines the boundary between different spin phases of the system ground state, 
decreases with increasing r. 
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II. MODEL AND METHOD 



The Hamiltonian for the coupled QD system is given by 

H = H orb + H z , (1) 

H orb = h( ri ) + h(r 2 )+C( ri ,r 2 ), (2) 

Mr) = ^(p + ^A) 2 + V(r), (3) 

C(n,r 2 ) = e 2 /e| ri -r 2 |, (4) 

H z = gf i B J2 B -Si- (5) 

i 

Here, we use the material parameters of GaAs, electron effective mass m* = 0.067m e , 
dielectric constant e = 12.9, and g-factor g = —0.44. /ig is the Bohr magneton, and A = 
|[— By, Bx, 0] is the vector potential for the constant magnetic field B oriented perpendicular 
to the QD plane (xy-pl&ne) . The Zeeman effect simply induces a lowering of the single- 
particle (SP) and triplet energies by 13 and 25 /xeV/T, respectively. 
We use the following model potential for the coupled QD system:- 

= _y Le -(x+d/2f/Rl+y 2 /Rl 

where Vl and Vr are the depth of the left and right QDs (equivalent to the QD gate voltages 
in experimental structures^) which can be independently varied, d is the interdot separation, 
R x and R y are the radius of the each QD in the x and y direction, respectively. In this work, 
we fix R x = 30 nm, and define QD aspect ratio r = R y /R x . Numerical exact diagonalization 
technique is used to solve for the single- and two-electron energies. Details of the method 
are published elsewhere.- 1 ^ 

Upon completion of the diagonalization procedure, we extract the SP energies and 
the two-particle energies E- /T . Here, "5" ( "T" ) denotes the singlet (triplet) state (In this 
paper, if not otherwise mentioned, "singlet" and "triplet" refer to the singlet and triplet 
states lowest in energy, respectively). The chemical potential of the N-th electron is given 
by the following equation: 1 

H(N) = E (N)-E (N-1), (7) 
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where Eq(N) [note E o (0) = 0] refers to the ground state energy with N electrons in the 
system. The exchange energy is given by 

J = £ T (2)-£ S (2). (8) 

For further analysis, the total energy of the two-electron system is partitioned into the 
expectation values of the SP energy K and Coulomb energy C 

E S ' T = (* S /T \H\* S /T ) 

= (^ /T |/ i (r 1 )+/ i (r 2 )|* S/T ) 

+ (^ /T |C( ri) r 2 )|^ /T ) 
= K S ' T + C 5 / T , (9) 

while the spectral function is defined as the projection coefficients of the lowest singlet and 
triplet states onto the SP product states^ 

«2 T = (^(ri)^(r 2 )|* 5/T (ri,r 2 )) . (10) 

The electron density is given by 

P 5/T (rx) = / |^ /T (r 1; r 2 )| 2 dr 2 . (11) 

Finally, the expectation value of the parity operator is given by 

^ S Q /T {-x u - yi ,-x 2 ,-y 2 )), (12) 

and for the parity operator along the |/-axis 

(pf T ) = (^ S /T {x 1 ,y 1 ,x 2 ,y 2 ) 

^o /T {xi,-yi,x 2 ,-y 2 )\ . (13) 

III. RESULTS 

A. Aspect ratio dependence of the exchange energy 

Figure [1] top panels show the potential contour plots r = 1 (left), r = 4 (middle), and 
r = 8 (right). As r increases, the potential becomes more elongated in the ^/-direction, while 
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the effective interdot distance (i.e., the x-distance between the two minima of the potential) 
and the interdot barrier height remain constant at 40 nm and 1.98 meV, respectively. 

In the lower panel of Fig. [TJ we plot the three lowest singlet (red/gray, solid) and triplet 
(blue/dark gray, dashed) energy levels as a function of r. With increasing r, the SP energies 
decreases (not shown), resulting in the decrease of the two-particle energy levels. We note 
that the lowest energy of the singlet state [Eq (2)] decreases smoothly with r, while the 
lowest energy of the triplet state [Eq (2)] exhibits a cusp at r = 3.9 because of the crossing 
of the lowest two triplet state energy levels. This cusp results in a sharp variation in the 
exchange energy dependence on r, which is shown in the inset of the lower panel of Fig. [U 
In the same inset, we show the variation of the tunnel coupling 2t = e\ — e§. For r < 4.3, 
the SP ground and first excited states have s and p x characters, respectively, and 2t barely 
increases from 1.8105 to 1.8114 meV with increasing r, because the energy contributions 
from the y-direction to eo and e\ cancel out. For r > 4.3, the SP first excited state bears a 
p y character, which causes 2t to decrease monotonically with r. 

In order to investigate in detail the cusp in the lowest triplet state energy, or, the crossing 
between the two lowest triple levels in the lower panel of Fig. [U we plot in Fig. [2] the spectral 
function of the two-electron wavefunction. It is seen that at r = 3.9 the triplet mainly 
consists of the [1, 2] and [2, 1] SP state pair, while at r = 4 it mainly consists of the [1, 3] 
and [3, 1] SP state pair. Here, 1, 2 and 3 denote the SP states in ascending energy, which 
have s, p x and p y characters, respectively, as shown in the Fig. [2] inset. Since the energy 
ordering of these SP states does not change as r changes from 3.9 to 4 (not shown here), 
the cusp in the lowest triplet state is due to a sudden transition of the triplet wavefunction 
from occupying an sp x pair to an sp y pair. 

In Fig. [31 we plot separately the (a) SP < K > and (b) Coulomb < C > contributions to 
the singlet and triplet state energies as a function of the QD aspect ratio r. As r increases, 
the general trend for all these energy terms is to decrease, leading to decreasing singlet 
and triplet energies shown in Fig. [TJ For the singlet state, both < K > and < C > terms 
decrease smoothly with r. For the triplet state, however, a discontinuity is seen from r = 3.9 
to 4: < K > (< C >) suddenly increases (decreases) by 0.128 (0.607) meV. It now becomes 
clear that the transition of the SP configuration shown in Fig. [2] from the sp x pair to the 
sp y pair is favored by the lowering of the Coulomb interaction despite the increase in the 
SP energy. The insets in Fig. E^a) and (b) show that the difference in the Coulomb energy 
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between singlet and triplet states (AC = (C T ) — (C 5 ) < 0) is always overcome by the SP 
energy contribution (AK = (K T ^ — (K s ^ > 0), leading to a positive exchange interaction 
(J = AK + AC, see Fig. [1]). 12 The comparison between AK and the tunnel coupling < It > 
in the inset of Fig. [3]^a) shows that the SP energy contribution to the singlet and triplet 
states is strongly influenced by the Coulomb interaction and is quite different from the 
noninteracting picture. 

As a consequence of the sudden change in the SP occupation, the y-symmetry P y of the 
two-electron wavefunction of the lowest triplet state changes abruptly from 1 to —1, which is 
validated by direct calculation of P y . We point out that the crossing between the lowest two 
triplet states by increasing r is allowed because they possess opposite ?/-symmetry, which 
exemplifies the general von Neumann- Wigner theorem relating the molecular energy levels 
to the two-electron wavefunction symmetry.— 

The contour plots in Fig. H] clearly show that from r = 3.9 (first row) to r = 4 (second 
row) the electron density in the lowest singlet state barely changes, while the density in 
the lowest triplet state changes abruptly from two peaks localized in the left and right QDs 
(the separation of two peaks in the x-direction is ~ 40 nm) to four peaks separated along 
both x and y directions (separation between peaks in the x and y directions are 20 and 
40 nm, respectively), again due to the sudden change in the SP configuration. The third 
row in Fig. [3] shows that at r = 8, both the singlet and triplet densities exhibit four peaks 
separated in both the x and y directions. Our analysis shows that from r = 4 to r = 8, the 
left and right peaks in the singlet state density gradually separate into four peaks, and the 
separation between the top two and bottom two peaks in the triplet state density smoothly 
increases. Such electron localization effects at large r are discussed for other many-electron 
QD systems with weak confinement, see, e.g. Ref. [l4j and references therein. 

B. Stability diagrams 

In Fig. [51 upper panels, we plot the stability diagrams 1 of the coupled QDs for r = 1 
(left), r = 3 (middle), and r = 5 (right) for R x = 30 nm, d = 50 nm and B = T. The 
solid curves indicated by arrows shows the computed contours, where chemical potentials of 
the first electron (red), the second electron in the singlet state (green), and second electron 
in the triplet state (blue) are equal to the reference value [/i(l) = At s (2) = /U T (2) = —21 
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meV]. According to the general shape of the stability diagram for coupled QDs,- we use 
dotted straight lines on the diagrams to separate different charge states indicated by discrete 
electron numbers on the left and right QDs, e.g., (0, 1) means zero electrons on the left QD 
and one electron on the right QD. Specifically, the boundaries between the (1, 1) and (0,2) 
[or (2, 0)] states are taken extending from the point on the \i (2) curve at which the curvature 
is the largest for Vl ^ Vr, e.g., point C on the upper left panel and parallel to the main 
diagonal. In the absence of mag netic field (B = 0), the fi s (2) curve is the boundary between 
one and two electrons in the system (in the linear transport regime wherein the source and 
drain chemical potentials are nearly the same).— Based on this fact, we extrapolate from 
the first off-diagonal triple point (e.g., point C on the upper left panel) to get the boundary 
between two- and three-electron states [green dotted curve indicated by ;u(3)]. Here, we 
assume that the triple point separation between charge states (1,0) and (2, 1) (or between 
the (0, 1) and (1, 2)) is the same as the separation between the (0, 0) and (1, 1) states.— 

In Fig. [51 we notice that, on the one hand, as r increases, the crossing points of the /x(l), 
H (2) and [i T (2) curves with the main diagonal shift to smaller Vl = Vr values because 
the SP energies decreases as r increases, and as such a less negative Vl and Vr value is 
required to charge the coupled QDs. On the other hand, the double-triple point (DTP) 
separation, i.e., the separation between the crossing points of /i(l) and /i(2) curves with 
Vl = Vr, decreases with r. For the singlet (triplet) state, the DTP separation measured in 
AV L = AV R is 5.181 (5.269), 4.128 (4.725) and 3.473 (3.907) for r = 1, 3 and 5, respectively. 
This decreasing trend of the DTP separation suggests that the coupling strength between 
the two QDs decreases with increasing r (see Ref. 111). However, from our direct calculations 
shown in the inset of Fig. [TJ lower panel, the tunnel coupling decreases only for r > 4.3, 
while the exchange energy is largest for r = 3.9. The discrepancies regarding the coupling 
strength between the DTP separation and direct calculations of the tunnel and exchange 
couplings can be understood by observing the following: the DTP separation is given by 
2t + C, where 2t and C denote tunnel coupling and interdot Coulomb interaction. As r 
increases 2t decreases for r > 4.3, while < C > monotonically decreases for both singlet 
and triplet (see Fig. [3]). As a result, the DTP separation decreases. The exchange coupling, 
however, is determined by the energy difference between the singlet and triplet states. As 
shown in Fig. [31 such energy difference, when splitted into the SP contribution < AK > 
and the Coulomb contribution < AC >, has a complicated dependence on r. In contrast, 
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if the interdot separation were increased to decouple the two QDs, then all quantities 2t, 
C, < AK >, and < AC > would decrease, leading to both decreasing DTP separation and 
exchange energy.-*^ 

One important feature shown in Fig. is that as r increases, the distance between the 
triple points on the main diagonal and the first off diagonal (e.g., points B and C in the 
upper left panel of Fig. [5]) becomes smaller, and at large r these triple points coincide. 
Consequently, the (1,1) stability region shrinks and finally disappears. This is because at 
large aspect ratios, even a small amount of interdot detuning can localized both electrons 
into the lower QD, resulting in an unstable (1, 1) charge state. The boundary /i T (2) at r = 5 
suggests that the (1, 1) charge state is also unstable for the triplet state, although the /i S (2) 
and /x T (2) curves evolve in different fashion as r increases. 

After locating the different charge stable regions on the stability diagram, we now inves- 
tigate the interdot detuning effect by departing from the center of the (1, 1) region along 
the direction perpendicular to the main diagonal, i.e., Vl + Vr — constant. Such detuning 
effects are important as two electrons transfer to a single QD, which is a key step in spin co- 
herent manipulation and spin-to-charge conversion in two-electron double QD experiments 
for quantum logic gate applications.- 1 ^ 

The solid curves in Fig. [5j lower panels, show the exchange energy J as a function of 
interdot detuning e = Vj, — Vr along the Vl + Vr = constant line [e = is chosen at the (1, 1) 
region center]. In the case of coupled circular QDs (r = 1), both singlet and triplet states 
localize progressively into the lower QD with increasing e, leading to a monotonic increase 
of J. Such a dependence is similar to recent experimental and theoretical results. For 
r = 3, a sharp cusp in J occurs at e ~ 4 meV before which J monotonically increases with e. 
This cusp is induced by a sudden SP configuration change in the lowest triplet state, which 
is similar to the effects seen in Fig. [1] and analyzed in Fig. [2], albeit here the perturbation 
in the Hamiltonian is introduced by interdot detuning instead of deformation effects. More 
detailed analysis of the two-particle energies and electron density for the r = 1 and r = 3 



cases can be found in Ref. 



17 



For r = 5, we observe that the exchange energy decreases 
monotonically with e, because the Coulomb energy difference between the singlet and triplet 
states becomes smaller as the two electrons in both the singlet and triplet states localize at 
the opposite ends of the lower single QD to minimize their Coulomb interaction. 

In the lower panels of Fig. we also plot the e dependence of A ST (dashed curves), the 
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difference between the /x 5 (2) and /i T (2) curves projected along the main diagonal. A ST is 
relevant in this context because in coupled QD experiments the chemical potential contour 
lines are mapped out by single-electron charging measurements, which provides useful infor- 
mation on the electronic structure of the QD.-^'iS Here, we notice that although the general 
detuning dependence is similar between J and A ST , a linear factor is not sufficient to scale 
values of J to overlap with those of A ST because the two quantities are extracted under 
different bias conditions. It should be pointed out that transport experiments measure the 
quantity A ST , which differs quantitatively from the exchange energy J. 

In Fig. El we plot the charge stability diagram of the coupled QDs for r = 1 (left), 
r = 3 (middle), and r = 5 (right) for R x = 30 nm, d = 60 nm and B = T. Compared 
to the data in Fig. \5[ which correspond to strongly coupled QDs, the data in Fig. [6] depict 
the situation in decoupled QDs.- In this case, as r increases, (1) the crossing points of the 
//(l), // s (2) and /i T (2) curves with the main diagonal shift to smaller Vl = Vr values; (2) 
the DTP separation decreases [The DTP separation is 2.847 (2.860), 2.618 (2.641), 2.514 
(2.545), for r = 1, r = 3, and r = 5, respectively]; and (3) the (1, 1) region becomes smaller. 
These behaviors are similar to those for d = 50 nm. However, the (1, 1) region does not 
vanish at d — 60 nm and r = 5 because as the QDs are more decoupled, both the interdot 
distance and interdot barrier height become larger, which require a larger interdot detuning 
to "push" both electrons into the lower QD. At a fixed r, the DTP separation (curvature 
at the triplet points) is smaller (larger) for d = 60 nm than for d = 50 nm, indicating that 
both tunnel coupling and Coulomb interaction are smaller for more decoupled QDs. 7 

C. Spin phase diagram 

In this subsection, we discuss the variation of the exchange energy J as a function of both 
r and B. By identifying the regions where J assumes different signs, we construct the spin 
phase diagram in which the two-electron ground state spin state (either 5* = or 5* = 1) is 
shown function of r and B.— ^ 

In Fig. [7J we plot the exchange energy J as a function of the QD aspect ratio r and 
the magnetic field B perpendicular to the xy-plane. At fixed r, as B increases, J decreases 
from its value at B = T to become negative and saturate at very large magnetic field, as 
previously reported. 6 - We note that at intermediate r (r ~ 4), J changes much faster with B 
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than at small or large r. This B-field effect at intermediate r values is associated with the 2D 
confinement of the QDs, i.e., near r = 4 the SP level separations in the x- and y-directions 
are comparable (cf. Fig. [TJ lower inset, 2t curve). We also note that, with increasing r, the 
relative change of J is small for B ~ 1 T, while it is much larger for B~0TorB~2 
T. The kink in J at B = T (cf. Fig. [TJ lower inset, J curve), due to the crossing of two 
lowest triplet levels, does not exist for B ^ T because a nonzero magnetic field couples 
the SP states with different Cartesian symmetries, thereby removing the condition for the 
crossing of the lowest two triplet states. In the investigated ranges of r and B, J assumes a 
maximum (minimum) value of 0.773 (—0.372) meV at r 3.9, B m T (r 4.4, B m 1.6 
T). 

The projected contour plots in Fig. [7] shows that the first singlet-triplet transition (at 
which J first crosses zero as B increases from zero at fixed r) occurs at a smaller B value 
as r increases, which is shown by the thick white dashed curve on the contour plot in Fig. 
[3 Such a dependence can be understood by observing that, in the absence of the B field, 
as r increases the SP energy spacing decreases, and, for a larger r, a smaller magnetic 
field is needed to further decrease the SP spacing and bring the triplet state to the ground 
state with the aid of the Coulomb energy difference between the singlet and triplet states. 
At higher magnetic field and larger r, we observe another contour line for J = (thick 
solid white curve at the lower left corner). The reappearance of the singlet state as the 
ground state is reminiscent of the singlet-triplet oscillation found for a two-electron single 
QD and also reported elsewhere for two-electron QDs with strong confinement.— 1 ^ 1 ^ In the 
foregoing discussion, we had not included the Zeeman energy for the triplet state, which 
would lower the triplet energy such that the boundary for the first singlet-triplet transition 
(thick white dashed curve) would shift to lower values of r and B, while the second singlet- 
triplet transition (thick white solid curve) would move to higher values of r and B. 

IV. CONCLUSIONS 

We have shown that the exchange energy between two electrons in coupled elongated 
quantum dots is enhanced by increasing the aspect ratio of the dots in the perpendicular 
direction to the coupling direction. However, there is an optimum aspect ratio beyond 
which the electron density in each dot starts to localize, and the exchange energy decreases. 
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With increasing aspect ratio, the (1,1) region becomes unstable with respect to interdot 
detuning, which is undesirable for two spin-qubit operations. We have also shown that the 
exchange energy in symmetrically biased coupled quantum dots is tunable between maximum 
(positive) and minimum (negative) values by varying the magnetic field and the QD aspect 
ratio. 
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Figures 




FIG. 1: (Color online) Top panels: potential contour plots of coupled QD with r = 1 (left), r = 4 
(middle), and r = 8 (right). Redder (darker gray) regions correspond to higher potential. Bottom 
panel: three lowest singlet (red/gray, solid lines) and triplet (blue/dark gray, dashed lines) energy 
levels as a function of QD aspect ratio r. The inset shows r dependence of the exchange energy 
J (blue/dark gray, solid) and tunnel coupling 2t (red/gray, dotted). For all panels, Vl = Vr = 25 
meV, d = 50 nm, B = 0T. 
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FIG. 2: (Color online) Spectral decomposition of the two-electron wavefunction onto different 
single-particle pairs. The red (dark) columns are for r = 3.9, while the yellow (bright) columns 
are for r = 4. The inset shows the contour plots of lowest three single-particle states in ascending 
order (indicated by number) of energy for both r = 3.9 and r = 4. The state symmetry is shown 
in parenthesis. 
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FIG. 3: (Color online) (a) Single-particle energy contribution (K) as a function of QD aspect ratio 
r. The red/gray, solid (blue/dark gray, dashed) line is for the singlet (triplet) state. Inset: the 
solid line shows the difference AK = (K T ^ — (K s y The dashed line shows 2t as a comparison, 
(b) Coulomb energy contribution (C) as a function of QD aspect ratio r. The red/gray, solid 
(blue/dark gray, dashed) line is for the singlet (triplet) state. Inset: the solid line shows the 
difference AC = (C T ) - (C s ). 
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FIG. 4: (Color online) Contour plots of the electron density for both singlet (left column) and 
triplet (right column). Rows I, II and III are for r = 3.9, r = 4 and r = 8, respectively. In the 
plots, redder (darker gray) regions correspond to lower electron density. 
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FIG. 5: (Color online) Top panels: stability diagrams for r = 1 (left), r = 3 (middle) and r = 5 
(right). In each diagram, the red, green and blue curves (solid) are computed contour lines at 
which the chemical potential //(l), \i (2) and /i^(2) equal to the reference value \i re f = -21 meV, 
respectively. Curves for different chemical potentials are also indicated by arrows. The dotted 
straight lines are a guide for eyes separating different stable charge states. Note that the exact 
locations of the yu(3) curve (green dotted curve) and (1, 2), (2, 1) regions are not computed. In the 
left two top panels, the (1, 1) region is indicated by the shaded area. In the left most upper panel, 
we also indicate the double-triple points A and B. Point C is where the \x (2) curve has the largest 
curvature for Vl ^ Vr. For corresponding QD aspect ratios, the bottom panels show J (solid 
curves) as a function of interdot detuning e = Vl — Vr from the center of the (1,1) region. The 
dashed curves on the bottom panels show the separation (A ST ) between the contour lines of /i S (2) 
and l^ T {2) projected along the main diagonal of as a function of interdot detuning e = Vl — Vr. 
All data are obtained at R x = 30 nm, d = 50 nm and B = T. 
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FIG. 6: (Color online) Stability diagrams for r = 1 (left), r = 3 (middle) and r = 5 (right). In each 
diagram, the red, green and blue curves (solid) are computed contour lines at which the chemical 
potential n s (2) and /U T (2) equal to the reference value fi re f = —19 meV, respectively. Curves 
for different chemical potentials are also indicated by arrows. The dotted straight lines are a guide 
for eyes separating different stable charge states. Note that the exact locations of the //(3) curve 
(green dotted curve) and (1, 2), (2, 1) regions are not computed. In each panel, the (1, 1) region is 
indicated by the shaded area. All data are obtained at R x = 30 nm, d = 60 nm, and B = T. 




FIG. 7: (Color online) Mesh (contour) plot of the exchange energy J as a function of QD aspect 
ratio r and the magnetic field B. The thick white curves (solid and dashed) on the contour plot 
correspond to J = 0. Total spin of the two-electron ground state is given in different regions. 
Redder (darker gray) regions correspond to lower J value. 
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